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Abstract 



It is investigated the behaviour of the "constants" G, c and A in the framework of a perfect 
fluid LRS Bianchi I cosmological model. It has been taken into account the effects of a c— variable 
into the curvature tensor. Two exact cosmological solutions are investigated, arriving to the 
conclusion that if q < (deceleration parameter) then G, c are growing functions on time t 
CO t while A is a negative decreasing function on time. 



^ '. 1 Introduction 
O 

In a recent paper (see pQ) we have investigated the behaviour of the "constants" G, c and A in the 
frame work of a bulk viscous LRS Bianchi I cosmological model where the effects of a c— variable 
t^j- | into the curvature tensor were taken into account. The main conclusion of such work is that G 

and c, under the physical restrictions (thermodynamics constrains), are growing functions on time 
while A is a negative decreasing function. In that paper, all the physical quantities depend of the 
viscous parameter and the thermodynamical restrictions bring us to determine all the parameters 
without any ambiguity, concluding that only the physical solutions are those which equation of 
state is u) = 1 (ultra stiff matter) and G and c rea growing functions on time while A behaves as a 
negative decreasing function. 

As we have showed in such paper for late time viscous, effects vanish in such a way that the 
model tends to a perfect fluid era. In this paper we try to study precisely this situation. Taking 
into account the effects of a c— variable into the curvature tensor we outline the field equations for 
a perfect fluid cosmological model with LRS Bianchi I symmetries. Under the assumed hypotheses 
we shall study two exact solutions. In this occasion we will not be able to determine the numerical 
values of the constants without any ambiguity in such a way that only imposing the condition q < 
(i.e. that the model accelerates) we arrive to the conclusion that G and c are growing while A is 
negative and decreasing. We believe that these solutions must be correct since in some way these 
solutions have to connect with the viscous one. 

The paper is organized as follows. In section 2 we outline the field equations and the assumed 
hypotheses in order to find the exact solutions. In section 3 we show the solution of the two studied 
models. 

In the first of them, the solution has been obtained under the hypotheses that the scale factor 
are follow the relationship X(t) ~ Y(t) n , this hypothesis is standard in this class of models. We 
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also assumed that the speed of light, c(t), follows a power law dependence on time t, i.e. c(t) = cot a , 
where a is a undetermined numerical constants. With these two hypotheses it is founded a singular 
exact solution to the field equations, but as we have mentioned above, in this case, we are only 
able to find poor restrictions to the numerical values. In order to see how each quantity behaves, 
we have plotted some cases, finding that the "constants" G and c can be growing as well as 
decreasing functions (with some restrictions) while A is a negative decreasing function on time t. 
In the same way we show that it is possible to find (fine tuning) cases which verify the relationship 
G/c 2 = const.. Only imposing some restrictions, like q < 0, it is founded that G and c must be 
growing while A behaves as a negative decreasing function. 

In the second studied solution founded under the hypotheses c(t) « Y(t) n and X(t) ~ Y(t) m , 
we shall show that it is non-singular. In the same way than in the first solution we shall find 
some restrictions to our model and we will plot some cases in order to see how the main quantities 
behave. As in the above solution we arrive to the conclusion that the most probable solution is 
those where G and c are growing functions while A is a decreasing function. 



2 The model and the hypotheses 

We consider a locally rotationally symmetric (LRS) Bianchi I (Hj-jlj) spacetime with metric 

ds 2 = -c 2 (t)dt 2 + X 2 (t)dx 2 + Y 2 {t) (dy 2 + dz 2 ) , (1) 
where the gravitational field equations with variable G, c and A are: 



1 87rG(t) 

Rik - -^QikR = c4 ,^ T ik + A{t)g ik , (2) 



where the energy momentum tensor is: 

T ij = (P + P) UiUj - pg ij} (3) 

and where p and p satisfy the usual equation of state, p = ojp in such a way that uj £ (—1,1] , that 
is to say, our universe is modeled by a perfect fluid. 

Applying the covariance divergence to the second member of equation (j2J) we get: 

div[^Ti + 8i^j =0, (4) 

hence 



that simplifies to: 



rpj GA • C^jAj 

T ^-\~~~G) li ~~8^G~> (5) 



Ac 4 G c 
P + ^ + p)H = -^G- p G + ^ (6) 



Therefore, with all these assumptions and taking into account the conservation principle, i.e., 
div(T-) = 0, the resulting field equations are as follows: 
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Since we have defined the 4— velocity as: 

1 



u 



c(t) 



,0,0,0 



(7) 
(8) 
(9) 

(10) 

(11) 

(12) 



then the expansion is defined as follows 



and the shear is 



a 





(13) 



(14) 



Curvature is described by the tensor field Rj kt - 



It is well know that if one uses the singular 
behaviour of the tensor components or its derivates as a criterion for singularities, one gets into 
trouble since the singular behaviour of the coordinates or the tetrad basis rather than the curvature 
tensor. To avoid this problem, one should examine the scalars formed out of the curvature. The 
invariants K\ and K 2 (the Kretschmann scalars) are very useful for the study of the singular 
behaviour: 



R ijkl = 

c 4 




YYc 



(15) 
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2.1 Simplifying hypotheses 

In order to solve these equations we shall need to make the following simplifying hypotheses: 
1. From the equations (jHJ) and Q it is obtained 

Y Y 2 f X y\ c XY X n 

y + ^ + {x-y)c-xy-x=°- (17) 

(a) Then we can make the next hypothesis: we impose that 

X « Y n , (18) 

therefore equation (|17j) yields 

\ 1 — n J Y c 

obtaining a trivial solution 

Y = (c 2 (n + 2) + C^n + 2) J c(t)dt\™ , (20) 

then we will need to make a hypothesis about the behaviour of c(t) in order to obtain a 
complete solution. 

(b) If we make the next assumption 

c(i) = c Y n , (21) 
then eq. (|17jl has the following solution: 



2 — n 



Y = X I C 2 (n - 2) y X n ~ 6 dt + Ci(2 - n) 1 , (22) 

then we will need to make a hypothesis about the behaviour of X(t) in order to obtain 
a complete solution. 

In the next section we shall study two different solutions which correspond to each one of the 
hypotheses. 



3 Solutions 



3.1 Model 1. Hypotheses la. 

If we follow the hypotheses la then it is observed that 

Y = (c 2 (n + 2) + d(n + 2) J c(t)dt\ ^ , (23) 

and making the following hypothesis about the behaviour of c(t) 

c(t) = c t a , (24) 

with u£K, which is possibly the most simple hypothesis about it behaviour. We have choice this 
form taken into account the results obtained in our previous paper (see j7j) where we have founded, 
through the Lie analysis of the differential equations, that the only physical solution for the function 
c(i) follows a power-law dependence of time t. We can choice other form for this function and with 
the help of eq. (j23j) it is obtained the exact form of the scale factors Y{t) and X(t), which brings 
us to a very different scenario. Therefore, for our election of c(t) (J2U we found that Y(t) and X{t) 
behaves as 

Y=(c 2 (n + 2) +Cl ^M^, =, X=(c 2 (n + 2) +Cl ^t^ 
V a + 1 J \ a+l 

(25) 

observing that we must to impose the following restrictions on the possible values of the constants 
a and n, in this way, we see that a ^ —1 and n > —2 if we want that our scale factor Y(t) be a 
growing function. 
Now defining 

fx Y\ ,Y (a + l)Cit a , , 

g -^ + 2 Fj=' 2 + "»F- (c 2 (l + i) + c 1 ^) - <26> 

we can calculate the deceleration parameter q 

d /3\ . Co _ a _! 1 



(n + 2) 



as we can see this expression only has sense if C 2 = 0, therefore 

2 - a 



a + l' 



(28) 



as it is observed we only have q < if a > 2. 
Now taking into account the expression 



p + (p + p) H = 0, => p = d ^ (n+2)( " +1) , (29) 
hence the energy density behaves as 

p = d Q t-^ +1 \ (30) 
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being do = const. 
Now from eq. 



x Y 2 8ttG , , 
(2n + l)^ = -^p + Ac 2 , (31) 



and denning 

(2n+l)^ = /(*), (32) 

we obtain A as 

A =(/W-^"W)^). < 33 > 

and substituting this expression into 

we end obtaining the behaviour of the "constant" G as: 

g = S(H)> (35> 

therefore it yields 

G = G t 3a+a;(a+1) -\ (36) 
being Go = const.. As it is observed, the quotient between G and c behaves as 

% « t "-M»+i)-i (37) 
cr 

finding that we only have the relationship 

G l-w , , 

= const. <=r- a = t r-. (38) 

cr (1 + a; J 

We insist into emphasize this relationship, because in ([7]) we obtain this expression as integration 
condition (with out any supplementary hypothesis) i.e. at least this expression could have math- 
ematical sense (as symmetry condition). In a subsequent paper (|T), where we were studying a 
LRS Bianchi I viscous model, this relationship was also verifies iff the viscous parameter was the 
usual one i.e. 7 = 1/2. In the next subsection (see below) this relationship help us to determine 
the values of the numerical constants. 

Once it has been obtained G then we back to eq. (|33|) obtaining in this way the behaviour of 

A : 

A = A ^ 2(a+1) , (39) 

where Ao = const. 

The behaviour of the Krestchmann scalars are: 

^ 1 „ t -4(«+i) j K 2 *r< a+1 \ (40) 
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while the expansion and the shear behaves as: 



^±l t -( a+1 \ (41) 
Co 



* = ^<pz£r a + l)t-M (42) 
3 (n + 2) 

It is observed that all the quantities depend only of the constants a and u and only X(t) depends 
of n. Therefore from these results we only can say that n > —2, uj > —1 and a > — 1, we have not 
more restrictions for these constants except that if we want that our solution verifies the condition 
q < then o > 2, but this condition tells us that c(t) is a growing function on time t. 



3.1.1 Conclusion for this model. Numerical values and graphics for the main quan- 
tities. 

In the viscous case discussed in ( ) we were able to find a concrete physical values for the constants 
(a,n,uj) under thermodynamical considerations but in this case we only have a poor restrictions, 
n > —2, uj > — 1 and a > — 1 but if q < then a > 2. But we cannot to discard the case a = —1/2, 
this possibility is so valid than the case a = 1/2. In this subsection we will study numerically some 
of the possible cases in order to show the behaviour of the different quantities. In the solutions 
exposed above we have written, for example Go or Ao as constants, but these numerical constants 
are horrible expression that depend of the constants (a,n,uj) and other numerical factors. The 
next figures are sensible to all these numerical factors. We shall study five cases, for this purpose 
we need to fix some numerical constants (a, n) as well as consider different equations of state (to). 
Since we only have the above restrictions (which do not help us) we can choose these numerical 
values in a very arbitrary way but we prefer to "assume" that eq. (|38|) is verified (since in previous 
works this relationship was obtained as integration condition and hence we believe that our model, 
and without other restrictions, must verifies it) in this way we find that the constant a = ip(uj) 
i.e. is a function of the equation of state and to fix arbitrarily the value of n since it only affects 
directly to the scale factor X(t). 

In the rest of this section we shall consider the following values for the numerical constants 

C x = l, C 2 = 0, c = l, do = l, 

and in the following table is summarized the colour of each solution which correspond to the different 
values of (a, n, uj) and the corresponding value of the parameter q : 



colour 


UJ 


a 


n 


Q 


red 
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-1/3 


3/2 


7/2 > 


blue 


1/3 


1/2 


1/2 


1 > 


magenta 





1 


1/5 


1/2 > 


black 


-1/3 


2 


5/2 





points 


-2/3 


5 


2/3 


-1/2 < 



We would like to emphasize that the red colour solution i.e. it has been calculated with the 
following numerical values (uj = l,a = — 1/3, n = 3/2) does not verifies the condition (|38(1 since 
for uj = l,a = and it has q > 0, the rest of the solutions verify such relationship. The only 
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solution that verifies the condition q < corresponds to the equation of state uj = —2/3 which is 
in agreement with the recent observations and theoretical models. 

In the first place we begin studying the behaviour of the radius of our model. With these 
numerical values we can see in fig. Q the different behaviour of our scale factors. In all the 
studied cases we can see that these scale factors are growing functions on time t. These solutions 
are singular as the Krestchmann invariants show us. 




t i 



Figure 1: In the figure (a) it is plotted the radius Y(t) while in figure (b) it is plotted the variation 
of the radius X(t). Both radius have a singular start and are growing functions on time. 

With regard to the energy density (see fig. 0) we can see that all the solutions have physical 
meaning since this quantity depends only of the equation of state. 

10tt| 1 




Figure 2: The figure shows the variation of the energy density p(t). 

The variation of the "constants" G and c as well as the relationship G/c 2 is shown in fig. (J3J). 
These figures show us that both "constants" are growing functions on time t except in the case 
(red line) a = —1/3 where G and c are decreasing functions on time t, and as we have pointed out 
above this solution does not verify the condition G/c 2 = const.. The rest of solutions verify the 
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relationship G/c 2 = const, since we have chosen the numerical values with such purpose. 




1 2 3 4 5 

t 



Figure 3: In the figure (a) we show the variation of "constant" G(t). In the figure (b) we show the 
variation of the "constant" c(t) and in the figure (c) it is plotted the relationship G/c 2 . 

With regard to the cosmological "constant" , see fig. Q , it is observed that all the solutions are 
decreasing but negative. Maybe this fact tells us that we may consider A(i) as a true ghost energy 
density. 

The expansion and the shear behave as follows, see fig. ®. As we can see all the models studied 
show a decreasing expansion and the only models that have a positive shear are the plotted with 
the red and black color lines. The shear and expansion scalars calculated above indicate that the 
Universe is shearing and expanding with time. The shear, which is a degree of anisotropy in the 
Universe, decreases monotonically with time. This indicates the fact that the initially anisotropic 
Universe gradually tends to an isotropic Universe at late time (present epoch) which is in agreement 
with the recent observations of a negligible amount of anisotropy present in the CMBR. 

Other solutions can be studied changing the numerical values of the constants as well as taking 
into account different equations of state. We have found at least a solution in argument with the 
recent observations (q < 0) (black point) which equations of state are u = —2/3 (domain walls) 
with G and c as growing functions on time t and verifying G/c 2 = const.. 
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Figure 4: The variation of the cosmological "constant" A(t). 




Figure 5: In figure (a) it is plotted the expansion 9(t), while in figure (b) it is plotted the shear 
a(t). 



3.2 Model 2. Hypotheses lb. 

If we follow the hypotheses lb then it is observed that if we make the next assumption 

c(t) = c Y n , 

then eq. (fTTj) has the following solution: 



Y = X ( C 2 (n - 2) / X n - 3 dt + Ci(2 - n 



i 

2 — n 



(43) 



(44) 



then we will need to make a hypothesis about the behaviour of X{t) in order to obtain a complete 
solution. Now if we assume that 

X = Y m , (45) 
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with m 7^ l,then 

Y = (Ci (m - n + 2) t + C 2 (m - n + 2))(™-«+ 2 ) , 
with the restriction n < m + 2. Defining 



X V y v ; y (m - n + 2) (Ci* + G 2 ) ' 

we can calculate the deceleration parameter q 

d ( 3 \ 1 n 3n 
Q=~n 77 -1 = 2- 



dt\HJ m + 2' 

finding that 

g < 3n > 2 (m + 2) . 
Now taking into account the expression 

p + (p + p) # = 0, p = d Y~( 2+m ^ +1 \ 

we find the behaviour of the energy density 

-(2+m)(^ + l) 

p = d (Ci (m - n + 2) i + G 2 (m - n + 2)) (—"+2) , 



with do = const. 
Now from 



1N Y 2 8vrG . 2 
(271 + 1)—^ = -^p + Ac 2 , 



Y 2 c 2 
where 

(2n + l)^ = /(<), 

we obtain A as 

8^G , A 1 



m-—P(t)) wr A, 

and substituting this expression into 

Ac 4 G_ c 
87TG/9 + G c ' 

we end obtaining the behaviour of the "constant" G as 

, . % 2n+4+m+^(n + 2) 

G = *l ( L - | = G - (Cl (m " " + 2) 1 + C2 (m " n + 2)) (m "" +2) 



8vrp \J c) (Git + G 2 ) 3 

with Go = const.. As we can see, in this occasion, the relationship between the quantities 
behaves as: 

4+m+qj(re + 2) 

G (Gi(m-n + 2)t + G 2 (m-n + 2)) 

c 2 ~ (Git + G 2 ) 3 ' 
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and it is observed that 

G CCTlst .^ m= 3n + .(n + 2)_ lj 
cr 2 



we are following the same argumentation as in the above model. 

Once it has been obtained G then we back to eq. (|54[) obtaining in this way the behaviour of 

A : 



A = A, 



o 



— 2n 2(1 — ra) + u(n — m) 

(Ci(m-n + 2)t + C 2 (m-n + 2))(™-™+ 2 ) (d (m - n + 2) t + C 2 (m - n + 2)) (m-»+a) 



(d* + C 2 ) 2 (d* + d) 



(59) 



with A = const. 

The behaviour of the Krestchmann scalars are: 



4(n + l) 4(n + l) 

^ w (Cii + C 2 ) _ ^+2 , K 2 ^{C 1 t + C 2 T^+5 , (60) 

while the expansion and the shear behave as: 

m + 2 

9 = C 1 (m + 2)(C 1 t + C 2 f W ^ , (61) 



/q _ m + 2 

a = ^-d (m - 1) (C x t + d)" 7 ™ ■ (62) 
In this case we have the following restriction: 

n < m + 2, from Y(t), 

3n > 2 (m + 2) , if <? < 0, 

n > —1, from .fQ, 

m > —2, from and a, 

which are very similar to the obtained in the above model. Now if we "assume" that the condition 
(|58j) must be verified then it is obtained a new restriction 

3n + u (n + 2) 

m = 1. 

2 

3.2.1 Conclusion for this model. Numerical values and graphics for the main quan- 
tities. 

As in the above model and following the same argument, in this subsection we will study numerically 
some of the possible cases in order to show the behaviour of the different quantities. We shall study 
five cases, for this purpose we need to fix some numerical constants (n, m) as well as consider 
different equations of state (oS). Since we only have the above restrictions (which do not help us) 
we can choose these numerical values in a very arbitrary way but we will "assume" that eq. ()58|) 
is verified. 

In the rest of this section we shall consider the following values for the numerical constants 

d = l, d = 0, co = l, d = l, 
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and in the following table is summarized the colour of each solution which corresponds to the 
different values of (m,n,u) and the corresponding value of the parameter q : 



colour 


UJ 
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red 
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1/4 


1/2 


17/10 > 


blue 


1/3 


3/5 


1/3 


43/35 > 


magenta 





10/9 


2/3 


3/4 > 


black 


-1/3 


5/2 


5/3 


-0.04 < 


points 


-2/3 


3 


11/6 


-8/23 < 



We would like to emphasize that the black colour (line) solution i.e. it has been calculated with 
the following numerical values {ui = — 1/3, re = 5/2, m = 5/3) does not verify the condition (|5<S[) . the 
rest of the solutions verify such relationship. The only solutions that verifies the condition q < 
correspond to the equation of state ui = —1/3 and ui = —2/3 which is in agreement with the recent 
observations. 

With these numerical values we can see in fig. © the different behaviour of our scale factors. 
In all the studied cases we can see that these scale factors are growing functions on time t. These 
solutions are non-singular as the Krestchmann invariants show us. 




2 4 6 8 10 2 4 6 8 10 

t t 



Figure 6: In the figure (a) it is plotted the radius Y{t) while in figure (b) it is plotted the variation 
of the radius X(t). Both radius have a nonsingular start and are growing functions on time. 

With regard to the energy density (see fig. ((7J)) we can see that all solutions have physical 
meaning, since they are decreasing functions on time t. 

The variation of the "constants" G and c as well as the relationship G/c 2 is shown in fig. ©. 
These pictures show us that both "constants" are growing functions on time t. We would like 
to emphasize that only for the case u> = —1/3 we have that it is not verified the relationship 
G/c 2 = const, as we have imposed. 

With regard to the cosmological "constant", see fig. Q, it is observed that all the solutions 
are decreasing but negative. In these cases and because our models are non-singular A — ► ±const 
when t — ► 0. 

The expansion and the shear behave as follows, see fig. (|10j) . As we can see all the models 
studied show a decreasing expansion and only the models that have a positive shear are plotted 
with the red and blue colours. 
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Figure 7: Figure (a) shows the variation of the energy density p(t). Figure (b) shows the variation 
of the bulk viscous pressure II (t). 




0.6- 
0.4 : 



0.2 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 
t 

Figure 8: In the figure (a) we show the variation of "constant" G(t). In the figure (b) we show the 
variation of the "constant" c(t) and in the figure (c) it is plotted the relationship G/c 2 . 
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Figure 9: The variation of the cosmological "constant" A(t). 
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Figure 10: In figure (a) it is plotted the expansion 0(t), while in figure (b) it is plotted the shear 
a(t). 

Other solutions can be studied changing the numerical values of the constants as well as taking 
into account different equations of state. We have found at least two non-singular physical solution 
(black line and black points) which equations of state are u = —1/3 (strings) and u = —2/3 
(domain walls). 

As we have been able to see it is very difficult to determine the behaviour of the constants 
since we cannot determine or restrict better the value of the numerical constants that define such 
behaviour. While in the viscous case we arrive to the conclusion that G and c only could behave 
as growing functions on time t, in these models we have founded that a decreasing behaviour is 
allowed as well as the growing one. Only if we insist in that our model verifies the condition q < 
we find that G and c must be growing functions while A is a negative decreasing function on time 
t. 

We have chosen our numerical values under the "assumption" that G/c 2 = const, (fine tuning) 
but of course this restriction is arbitrary although in the viscous model was verified (only for 
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7 = 1/2, fine tuning again). With out other physical considerations we cannot ensure that G and 
c are growing or decreasing functions on time, nevertheless, it would seem reasonable to choose 
growing functions in order to connect the viscous era to the perfect fluid era. In both cases A is a 
negative decreasing function. 
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